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Abstract: We study the effect of magnetic field B on the critical temperature T c of 
the confinement-deconfinement phase transition in hard-wall AdS/QCD, and holographic 
duals of flavored and unflavored J\f = 4 super-Yang Mills theories on M 3 x S 1 . For all 
of the holographic models, we find that T C (B ) decreases with increasing magnetic field 
B <C T 2 , consistent with the inverse magnetic catalysis recently observed in lattice QCD 
for B < 1 GeV 2 . We also predict that, for large magnetic field B T 2 , the critical 
temperature T C (B), eventually, starts to increase with increasing magnetic field B T 2 
and asymptotes to a constant value. 
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1 Introduction 

Recently, the study of the QCD phase diagram for magnetic field B has attracted consid¬ 
erable attention [1-19], see [20] for a review. The main motivation for studying the QCD 
phase diagram under external magnetic field B stems from the fact that strong magnetic 
field B is produced in heavy ion collisions experiments at RHIC eB 0.01 GeV 2 and LHC 
eB ~ 0.25 GeV 2 [21], due to the charged spectator particles, which has interesting effects 
on the quark-gluon plasma created during these heavy ion collision experiments [22-27], 
see [20] for a review. A strong magnetic field eB ~ 4 GeV 2 is also produced during the 
electroweak phase transition of the early Universe [28], and relatively weaker magnetic field 
eB ~ 1 AleV 2 is produced in the interior of dense neutron stars [29]. 

Another motivation comes from the fact that the study of the QCD phase diagram 
with magnetic field B is amenable to numerical simulations of QCD on the lattice, without 
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facing the sign problem of lattice QCD that exist in the case of non-zero baryon chemical 
potential [i b , creating an opportunity to compare the holographic and low energy effective 
models of QCD directly with QCD itself. 

Regarding the study of the QCD phase diagram for magnetic field B , most of the mod¬ 
els for QCD [1-10], including the holographic ones [11—16], have studied chiral-symmetry- 
restoration transition and have predicted that the critical temperature T c of the transition 
increases with increasing magnetic field B at zero chemical potential fi = 0. This enhanc¬ 
ing effect of the magnetic field B on the critical temperature T c has been termed magnetic 
catalysis. However, recent lattice QCD result [19] has indicated the opposite effect, that 
is, the critical temperature T c decreases with increasing magnetic field B , for B < 1 GeV 2 
and zero chemical potential /x = 0. This inhibiting effect of the magnetic field B on the 
critical temperature T c has been termed inverse magnetic catalysis. 

Even though, the recent lattice QCD result [19] has also indicated that the confinement- 
deconfinement and chiral symmetry breaking phase transitions occur at the same critical 
temperature T C (B) at least for B < 1 GeV 2 , most holographic calculations so far [11-17] 
have been concerned only with T C (B) of the chiral symmetry breaking phase transition. 

However, recently, reference [18], inspired by the recent lattice QCD result [19], has a 
priori assumed confinement and chiral symmetry breaking transitions to occur at the same 
critical temperature T c in Sakai-Sugimoto model, and has argued that, in this case, T C (B) 
must be a decreasing function of B , consistent with the recent lattice QCD result [19], but 
has not provided a direct computation of T C (B). 

In this paper, we give a direct computation of the critical temperature T C (B) of the 
confinement-deconfinement phase transition in hard-wall AdS/QCD, and holographic duals 
of flavored and unflavored Af = 4 SYM on M 3 x S 1 where S 1 is a circle of length l in one of 
the spatial directions. (Note that, at finite temperature T, M 3 is really S 4 x M 2 where S'/ is 
the thermal circle with length /,.) Also, note that, since the fermions of both the flavored 
and unflavored Af = 4 SYM on M 3 x S 1 obey antiperiodic boundary conditions around 
the circle S 1 , they acquire a tree-level mass m ~ j. The scalars are periodic around the 
circle, hence they acquire masses only at the quantum level through their couplings to the 
fermions [34]. The gluons, however, do not acquire masses, therefore, at low-energy, both 
flavored and unflavored Af = 4 SYM on M 3 x S 1 reduce to pure 3D Yang-Mills theory. 

It is well known that both flavored and unflavored AT = 4 super-Yang Mills theories 
(SYM) on flat spacetime M 4 are not confining gauge theories. However, they can be made 
confining in the large- N c limit by placing them on a compact space with length l, and 
the confinement-deconfinement phase transition occurs at critical temperature T c = ) 
[30, 31, 33], see [34, 35] for a review. In our case, the compact space is M 3 x S 1 , that is, we 
compactify one of the spatial dimensions into a circle of length l. 

The confinement-deconfinement phase transition both in flavored and unflavored Af = 
4 SYM on M 3 x S 1 is holographically modeled by a phase transition between a black hole 
solution with radius of horizon r = rh, and AdSs-soliton solution which smoothly ends 
at r = tq. However, to study the confinement-deconfinement phase transition in QCD 
on M 4 at strong coupling, we use the hard-wall AdS/QCD model where the confinement- 
deconfinement phase transition, of QCD on M , is holographically modeled by a phase 


- 2 - 


transition between a black hole solution with radius of horizon r = r^, and thermal- Ad,S§ 
solution with hard-wall IR cut-off r = r$. 

We derive the corresponding thermal- AdS§ solution which is the holographic dual to 
the confined phase of QCD on M 4 by starting from a black hole solution, which corresponds 
to the deconfined phase of strongly coupled QCD on M 4 , by setting the mass of the black 
hole to zero [47]. And, we derive the corresponding AdSs-soliton solution, which is the 
holographic dual to the confined phase of flavored and unflavored Af = 4 SYM on R 3 x S 1 , 
by ’’double Wick rotating” a black hole solution [34, 35]. 

In this paper, we use two black hole solutions in the presence of constant magnetic 
held B. First, we use the black hole solution in the presence of constant magnetic held 
B <T“ found in [36] to study the conhnement-deconhnement phase transition in strongly 
coupled QCD on M 4 and unflavored Af = 4 SYM on l 3 x S 1 . Then, we use the black hole 
solution in the presence of constant magnetic held B, including the backreaction of Nf 
havor or D7-branes for Nf -C N c , found in [38] to study the conhnement-deconhnement 
phase transition in flavored Af = 4 SYM on l 3 x S 1 . 

The effect of magnetic held B on different observables has also been studied in [39-43] 
using the backreacted black hole solution of [36] without havor D7-branes. 

Depending on the specihc holographic models to QCD, various length and energy scales 
appear throughout this paper. Some of the relevant length and energy scales are: the radius 
of the AdS ;5 spacetime L which we set to L = 1, the radius of the black hole horizon r/,. 
which is related to the Hawking temperature Tu of the black hole (which is dual to the 
held theory temperature T = Th), the radial position of the canonical singularity of the 
AdSs-soliton ro = ttT c (B = 0) = 7 r x 0.175 GeV = 0.55 GeV for flavored and unflavored 
Af = 4 SYM on M 3 x S 1 , the radial position of the hard-wall ro = ^05 = 0-323 GeV in the 
thermal- AdS§ solution for the hard-wall AdS/QCD, and an external magnetic held B in 
the range of 0 — 0.35 GeV 2 for the hard-wall AdS/QCD model and 0 — 4.2 GeV 2 for the 
havored Af = 4 SYM on R 3 x S 1 . 

The outline of this paper is as follows: In section 2, we write down the 5-dimensional 
Einstein-Maxwell action (which will be used to study conhnement-deconhnement phase 
transition in hard-wall AdS/QCD and holographic dual of unhavored Af = 4 SYM on 
R 3 x S 1 ) including the Gibbons-Hawking surface term and the appropriate counter terms. 
We then review the black hole solution in the presence of constant magnetic held BcT 2 
found in [36]. Then, starting from the black hole solution, which corresponds to the 
deconfined phase of strongly coupled QCD on hat spacetime and unhavored Af = 4 SYM 
on R 3 x S 1 , we derive the corresponding thermal-AdS '5 and AdS' 5 -soliton solutions, which 
correspond to the confined phases of strongly coupled QCD on hat spacetime and unhavored 
Af = 4 SYM on M 3 x S 1 , respectively. We also determine the on-shell Euclidean actions 
(free energies) for the black hole, thermal- AdS §, and AdS's-soliton solutions. 

In section 3, we compute the critical temperature T c of the conhnement-deconhnement 
phase transition in hard-wall AdS/QCD, and holographic duals of havored and unhavored 
Af = 4 SYM on R 3 x S 1 . We hrst compute the critical temperature T c of the conhnement- 
deconhnement phase transition in hard-wall AdS/QCD by requiring the difference between 
the black hole and thermal- AdS§ on-shell Euclidean actions vanish at T = T c . Then, we 
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compute the critical temperature T c of the confinement-deconfinement phase transition in 
holographic dual of unflavored M = 4 SYM on M 3 x S 1 by requiring the difference between 
the black hole and AdSs-soliton on-shell Euclidean actions vanish at T = T c . Finally, 
using the insight we gained, in computing the critical temperature T c of the confinement- 
deconfinement phase transition in holographic dual of unflavored M = 4 SYM on R 3 x S 1 , we 
compute the critical temperature T c of the confinement-deconfinement phase transition in 
holographic dual of flavored M = 4 SYM on K 3 x S 1 by constructing the backreacted AdS 5 - 
soliton solution from the backreacted black hole metric of D3/D7 model, with magnetic 
field B, found in [38]. 

In Appendix A, we compute the critical temperature T C (B) of the confinement-deconfinement 
phase transition in hard-wall AdS/QCD for large magnetic field B T 2 . 


2 Einstein-Maxwell theory in 5D 


In this section, we review elements of Einstein-Maxwell theory in 5D which will, sub¬ 
sequently, be used to study confinement-deconfinement phase transitions in hard-wall 
AdS/QCD and holographic dual of unflavored M = 4 SYM on K 3 x S 1 . 

The action of five-dimensional Einstein-Maxwell theory with a negative cosmological 
constant is [36]* 

S = Sbulk + 'S’bndy ) (2-1) 

where the bulk action S/mik is 


Sbulk 


1 


167tG5 

and the boundary action S'bndy is 

1 


r d 5 xV=g(R ~ F mn F MN + 


Sbndy — 


87rG5 




( 2 . 2 ) 


(2.3) 


r=r A 


In the boundary action S'bndy (2-3), the first term is the Gibbons-Hawking surface term, and 
the other terms are the counter terms needed to cancel the UV(rA —> 00 ) divergences in the 
bulk action in accordance with the holographic renormalization procedure [46]. Note that 
the counter terms are entirely constructed from the induced metric 7 ^ on the boundary 
surface at r = r\, that is, 


y) = diag (g tt (r A )^xxir^^yyir a), g zz (r a)) . (2.4) 

And, K is the trace, with respect to 7 ^, of the extrinsic curvature of the boundary given 
by K lw = (d r -y liu )/(2^g^). Using the matrix formula d fl (detM) = detMtv(M~ 1 d f ^M) 
[35], we can write K = ■j IJ ' , 'K llL , = [35, 36]. 

In addition to the Bianchi identity, the field equations are [36] 

Rmn = -j^Smn ~ FpQg M N + 2 F M pF N p , (2.5) 

*Our conventions here are such that the Ricci scalar here Rhere is related to the Ricci scalar there Rthere 
(given in [36]) by R here = -Rthere- 
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( 2 . 6 ) 


V m F M tv = 0 . 

From now on we set the AdS radius to unity, that is, L = 1. 

Turning on a constant bulk magnetic field, in the ^-direction, B z = F xy = d x A y — 
d y A x = B, where the bulk gauge potential A y (x,r) = ^B(x5]{ — ySfyJ which solves 
Maxwell’s equation (2.6), and contracting Einstein’s field equation (2.5), one can find the 
Ricci scalar R = g AIN Rmn to be 

R = -20+-B 2 g xx g yy . (2.7) 

3 

So, the on-shell Euclidean action Se (which can be found from the Lorentzian action (2.1) 
by analytic continuation in the imaginary time direction, i.e., f_E = it) takes the form 

Se = S^uik + ‘S'bndy ) ( 2 - 8 ) 

where the on-shell Euclidean bulk action S^ ulk is 

SU = 8)^5 Jo dtE f^( 4 + \ b2 9 XX 9 VV ) , (2.9) 

and, the on-shell Euclidean boundary action Sb ndy is 

S’Liy = - 3 + B‘ 2 g xx gyy\nr lX ] , (2.10) 

and, r\ is the UV cut-off while r' is the radius of the horizon r' = r^ for a black hole 
solution, and IR cut-off r' = ro for a thermal- AdS§ or AdS's-soliton solutions. From now 
on we set V 3 = 87^5 = 1. Also, note that the on-shell Euclidean action Se is related to 
the free energy F by Se = f3F. 


2.1 Background solutions with B«T 2 

In this subsection, we review the black hole solution in the presence of constant magnetic 
field B<r ] found in [36] which corresponds to the deconfined phase of strongly coupled 
QCD on M 4 (flat spacetime) and unflavored AT = 4 SYM on M 3 x S 1 . Then, starting 
from the black hole solution, by setting the mass of the black hole to zero [47], we derive 
the corresponding thermal- AdS^ solution which is the holographic dual to the confined 
phase of strongly coupled QCD on flat spacetime M 4 . And, by ’’double Wick rotating” the 
black hole solution [34, 35], we derive the corresponding AdSs-soliton solution which is the 
holographic dual to the confined phase of unflavored and strongly coupled AT = 4 SYM on 
M 3 x S 1 . 

^Note that the bulk gauge potential A^(x,r) and the corresponding bulk magnetic field B = F xy (x,r) 
are dual to the boundary gauge potential A y (x) = A y ,(x, r = 00 ) and the corresponding boundary magnetic 
field B = F xy {x,r = 00 ) which couple to the (7(1) charged particles of the field theory living at the 
boundary. Later on, when we start discussing specific holographic models to QCD, we will specify the type 
of global (7(1) gauge group and the associated boundary current J^(x). 
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2.1.1 Black hole 


For BcT 2 and electric charge density p , the perturbative black hole solution in powers 
of B, up to an integration constant 03 is given in Eq. 6.1 and 6.16 of Ref. [36]. Here, 
we set the electric charge density p = 0 and fix the integration constant 03 = — | so that 
the perturbative solution in powers of B matches the near boundary solution which is also 
given in Eq. 4.4, 4.5 and 6.16 of [36]. Therefore, the black hole solution in Eq. 6.1 and 
6.16 of Ref. [36], for vanishing electric charge density p = 0 and 03 = — |, takes the form 


f(r) 

q(r) 

h{r) 


r 2 (- f{r)dt 2 + q(r)dz 2 + h(r ) (dx 2 + dy 2 )) + 2 , 


dr 2 

f{r)r 2 


^l^ + 0(B 4 ), 

i + \b 2] ^ + o^), 


( 2 . 11 ) 


and, the Hawking temperature T becomes 


T=l = U'(r h ) 


2vr 






+ 0(B 4 ) , 


( 2 . 12 ) 


where M is the mass of the black hole, U(r ) = r 2 f(r ), the radius of the horizon r^ is 
defined by requiring f(r = r^) = 0, T is the Hawking temperature of the black hole, and ft 
is the length of the thermal circle which acquired a fixed value as a function of in order 
to avoid the canonical singularity at the horizon r = r^. One can also check that (2.11) 
indeed satisfies the Einstein field equation (2.5) or its contracted version (2.7). 


2.1.2 Thermal-AdS ,5 

The thermal-HdS '5 solution can be found from a black hole solution by setting the mass of 
the black hole M to zero, see [47] for the electrically charged black hole case. Therefore, 
from the black hole solution for B<T 2 ( 2 . 11 ), we can determine the thermal-AdS solution 
for B <C A. 2 r ~ Tq by setting the mass of the black hole M = 0, 

dr 2 

^thermal = ^ ( -fo(r)dt 2 + q{r)dz 2 + h{r) ( dx 2 + dy 2 )) + » ( 2 - 13 ) 

AM = i - \ b + e>(B 4 ). 

A(r) = 1 + + OfB 1 ) . 

2.1.3 AdSs-soliton 

The AdS' 5 -soliton solution [32, 33] can be determined from the black hole solution (2.11) 
by ’’double Wick rotation” t = iz' and z = it' [34, 35]. Therefore, for B <C Aj R ~ Tq the 
AdSs-soliton solution is, 
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where l is the length of the circle in the compactified z' direction which is arbitrary for the 
black hole solution but in order to avoid the canonical singularity at r = ro (where ro is 


defined by requiring / s (r = ro) = 0 ), it acquires a finite value which is given in terms of ro 
for the ylcZiS^-soliton solution (2.14). 

2.2 On-shell Euclidean actions with 6<CT 2 

In this subsection, we determine the on-shell Euclidean actions (free energies) for the black 
hole, thermal-Olds's, and AdSs-soliton solutions. And, we compute the difference between 
the on-shell Euclidean actions of the deconfining geometry (which is the black hole geometry 
for both hard-wall AdS/QCD and holographic dual of unflavored Af = 4 SYM on M 3 x S 1 ) 
and the confining geometry (which is the thermal-AdS's geometry for hard-wall AdS/QCD, 
and the AdS' 5 -soliton geometry for holographic dual of unflavored Af = 4 SYM on M 3 x S 1 ). 

2.2.1 Black hole 

The on-shell Euclidean action Se = •S'bh (2-8) for the black hole solution with B <C T 2 


( 2 . 11 ) is 


(2.15) 



where the on-shell Euclidean bulk action of the black hole J S' bulk for B <C T 2 is 



(2.16) 


and the on-shell Euclidean boundary action of the black hole S'bndy for B<T 2 is 



(2.17) 


The bulk action S^uik (2-16) (after using the black hole metric for B « f 2 (2.11), using 


the fact that h(r)y/q(r) = 1 + 0(B 4 ), evaluating the integrals, and simplifying) become 



(2.18) 


- 7 - 







(2.19) 


which diverges when r A —> oo, and the boundary action S bn dy (2-17) becomes 
Sbndy = -Piri + \b 2 Inr A - \m - l -B 2 ) + 0(B 4 ), 

where we ignored terms which goes to zero in the rA —> oo limit. Also note that (2.19) 
diverges when r A —> oo, but the sum of S'buik (2.18) and Sbndy (2-19) is finite. Hence, the 
black hole on-shell Euclidean action Sbh (2.15) is 

Sbh = Sbnik + S bndy = -p (V 4 - l -M + 2 -B 2 In r h - *H 2 ) + 0(B 4 ) . (2.20) 

2.2.2 Thermal- AdS§ 

The on-shell Euclidean action Se = Sthermai (2-8) for the thermal-AdSs solution with 
B <C A‘j R ~ rg (2.13) is 


Sthermai — Stbulk ~k Stbndy) 


( 2 . 21 ) 


where the on-shell Euclidean bulk action Stbulk of the thermal- AdS$ for B <C A 2 R ~ r§ is 

rP' r r a / 9 \ 

Stbulk = J o dt E J dr ^(4 + ^B 2 g xx g yy j , (2.22) 

and the on-shell Euclidean boundary action Stbndy of the thermal- AdS§ for B <C T 2 is 

Stbndy = ~j* - 3 + B 2 g xx g yv In r A ) . (2.23) 

The thermal-yldSs bulk action S t b u ik ( 2 . 22 ) (after using the thermal-AdSs metric for 
B <C T 2 (2.13), using the fact that h(r)^/q(r) = 1 + 0(B 4 ), evaluating the integrals, 
and simplifying) becomes 

Stbulk = -/3'(r 4 -rl-^ 2 lnr A + ^ 2 lnr 0 ) + 0 (H 4 ), (2.24) 

which diverges when r\ —oo, and the thermal- AdS§ boundary action 5 t bndy (2.23) be¬ 
comes 

Stbndy = ~P '(4 + ^ 2 lnr A - in 2 ) + 0(B 4 ), (2.25) 

which diverges as well when r A —> oo. But, the sum of S tbu ik (2.24) and S t b n dy (2.25) is 
finite. Hence, the thermal on-shell Euclidean action Sthermai (2.21) becomes 

Sthermai = ~P{r 4 + \b 2 lnr 0 - ^B 2 ) + 0(B 4 ). (2.26) 

where we used /3' = (3y/~f = /3 at the boundary r A —> oo. 

Therefore, ASe (which is the difference between the AdS§ black hole ( 2 . 20 ) and 
thermal- AdS§ (2.26) on-shell Euclidean actions) becomes 

ASe = Sbh - Sthermai = ~P (r 4 h - r 4 - -M + 2 -B 2 ln(^)) + 0(B 4 ) . (2.27) 
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2.2.3 AdS, 5 -soliton 

Since, black hole (2.11) and AdS' 5 -soliton (2.14) are equivalent Euclidean geometries, their 
on-shell Euclidean actions take similar form. In fact, the on-shell Euclidean action of AdS§- 
soliton can be found by merely replacing ry by ro in the on-shell Euclidean action for the 
black hole [35]. Therefore, the difference between the on-shell actions S'bh of the black hole 
(2.20) and S^iiton of AdS^-soliton geometries is simply 

AS e = S bh - S soliton = -fi (r 4 h - r 4 + ^B 2 In ^ + 0(B 4 ) . (2.28) 

3 Confinement-deconfinement phase transition in holographic models of 
QCD for B < T 2 

In this section, we compute the critical temperature T c of the confinement-deconfiuement 
phase transition in hard-wall AdS/QCD, and holographic duals of flavored and unflavored 
M = 4 SYM on M 3 x S 1 . 

3.1 Confinement-deconfinement phase transition in hard-wall AdS/QCD 

For hard-wall AdS/QCD [36] t, we determine the critical temperature T C (B ) of the confinement- 
deconfinement phase transition by first determining the critical radius of the horizon 
r h = i"hc from the condition that the difference between the Euclidean actions for the 
black hole and thermal- AdS$ solutions vanish at rh = r^c, i.e., ASe(Di = rh c ) = 0. For 
B <C T 2 , requiring ASe^/ic) = 0 in (2.27), we find the constraint equation for the critical 
radius of the horizon rh c to be 

r hc + \b 2 ln(—) - 2 tq + 0(B 4 ) = 0, (3.1) 

which can be solved numerically for rh c {B,ro). Note that, we have fixed M = 2?’q in 
(2.27), so that (3.1) reduces to the constraint equation found in [48, 49] at B = 0, which is 
rL = 2 fq. Once we find the solution for r^ c from the constraint equation (3.1), we can use 
(2.12) to find T c = T(r^ = rhc , M = 2?/)). The plot of the numerical solution for T C (B, ro) 
for B <C T 2 is given in Fig. 1, and the numerical plot clearly shows that T C (B) decreases 
with increasing B <C T 2 in agreement with the inverse magnetic catalysis recently found 
in lattice QCD for B < 1 GeV 2 [19]. 

3.2 Confinement-deconfinement phase transition in holographic dual of unfla¬ 
vored J\f = 4 SYM on M 3 x S 1 

For the holographic dual of unflavored AT = 4 SYM on M 3 x S ls , we study the confinement- 
deconfinement phase transition by using the same Einstein-Maxwell action in 5 D as we used 

^For the hard-wall AdS/QCD model, the bulk magnetic field B = F xy (x,r) and the corresponding bulk 
gauge potential A M (a;, r) = | B(xS / — y 8 are dual to the boundary magnetic field B = F xy (x, r = 00) and 
the corresponding boundary gauge potential A ll (x) = A fl (x,r = 00) of the U(l)v subgroup of the SU(Nf) 
global flavor group of QCD. And, the boundary vector gauge potential A ll (x) couples to the boundary 
conserved vector current Jy(x). 

§ For the unflavored Af = 4 SYM case, the bulk magnetic field B = F xy (x, r ) and the corresponding bulk 
gauge potential A / 1 (x,r) = \B(x 8 y ^ — yS^) are dual to the boundary magnetic field B = F xy (x,r = 00) 
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Figure 1. Critical temperature T C (B) of the hard-wall AdS/QCD with r 0 = ^105 = 0-323 GeV 
[48]. Note: B = \/3B is the physical magnetic held at the boundary [50]. 


for the hard-wall AdS/QCD, and the analysis will be similar to the hard-wall AdS/QCD 
case but, for the unflavored Af = 4 SYM on M 3 x S 1 case, we compactify the black hole 
solution in the z direction into a circle of length l, and compare its free energy with the free 
energy of AdSs-soliton solution (2.14) instead of the thermal- AdS^ solution (2.13) that we 
used for the hard-wall AdS/QCD. 

It is easy to see from (2.28) that the critical radius of the horizon rh = r^ c at which 
AS'e/'c^ = rhc) = 0 is given by = rhc = A)- Therefore, using (2.12), the critical 
temperature T c = T(rh = rhc = tq) becomes, 






In r 0 \ 
) 


+ 0{B 4 ) 


1 

7 ' 


(3.2) 


Fixing M = r q so that we reproduce the correct B = 0 result T C (B = 0) = ^r, and fixing 
ro from the value of T c at B = 0, which we denote as T°, we can write (3.2) in terms of 
T c = ? as 

T c = T c °fl-(|-) j +0(B a ) (3.3) 


where we dehned the critical magnetic field B c = an( ^ ^ ra dius of the 

AdS spacetime. From (3.3), it is easy to see that T c is a decreasing function with increasing 
B CT 2 in qualitative agreement with the recent lattice QCD result [19]. 


3.3 Confinement-deconfinement phase transition in holographic dual of fla¬ 
vored J\f = 4 SYM on M 3 x S 1 

In the previous subsection, we have studied the confinement-deconfinement phase transition 
in the holographic dual of unflavored Af = 4 SYM on M 3 x S 1 using the backreacted black 
hole and AdSs-soliton geometries, from which, we can infer a simple prescription of finding 
T c in any backreacted black hole and AdSs-soliton based models. 

and the boundary gauge potential A^fx) = A M (a;,r = oo) of the (7(1) subgroup of the SU(4)r global 
R-symmetry group of Af = 4 SYM which couples to the boundary R-current J^. 
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The prescription is, first find the backreacted metric and the Hawking temperature 
T{rh ) of the black hole, then the critical temperature T c is simply given by T c = T{rh = ro) 
where ro can be fixed by the value of Tj? = T C (B = 0). 

Therefore, using this prescription, we can determine T c (Nf , B ) of the holographic dual 
of flavored M = 4 SYM on M 3 x S 1 . To this end we will use the backreacted metric of 
D3/D7 model given in [38] where the authors have also found the Hawking temperature 
B) including the backreaction of Nf Z?7-branes and magnetic field B^ to be, see Eq. 
3.1 o'f [38] II, 

r = 7 ( 1 + i^f ( 1 ’ 2 \/^f)) +0(W/JV ' )2) ' (3 ' 4) 

Since the on-shell Euclidean action of the black hole solution (including the backreaction 
of Nf 177-branes and magnetic field B ) has also been given in Eq. 3.14 of [38], in order to 
find the corresponding Euclidean action of the HdSs-soliton, all we need to do is replace rh 
by ro in Eq. 3.14 of [38]. Hence, the difference between the two on-shell Euclidean actions 
vanishes at the critical radius of the horizon rh = The = r o- And, using rh = The = ? ’o in 
(3.4), the critical temperature T c = T(rh = Th c = fo) of the confinement-deconfinement 
phase transition in flavored N = 4 SYM on M 3 x S 1 becomes 

= 7 f 1 + ( 4 - !)) (3.5) 

which can be written in terms of T® = T c (Nf = 0, B = 0) = — as 

f 1 - 2 f+hS)) + ° m/N * >2} ’ (3 - 6) 

where A h is the value of the ’t Hooft coupling fixed at the horizon rh , that is, A h = 
ATrg s e^ h N c where g s is the string coupling constant and </>(r) is the dilaton scalar field. 

Note that, for B = 0, (3.6) reduces to 

7 = T c (i - Nj + 0((Nf /N c ) 2 ), (3,7) 

which is in a qualitative agreement with the hard-wall AdS/QCD [51], functional renormal¬ 
ization group study of QCD [52], and lattice QCD [53] results which show that T c decreases 

^For the D3/D7 model, we use a bulk magnetic field which corresponds to the Kalb-Ramond two form 
field B xy (x,r). And, also note that for the DBI action of the probe D7 brane, the gauge invariant and 
physically significant field strength is given by TYn = F m n + B m n = d m A n {x) — d„Am(x) which is the 
sum of Maxwell’s field strength F mn and the Kalb-Ramond two form field B mn . And, the bulk magnetic 
field B xy (x, r) and the corresponding bulk gauge potential Ay(x, r) are dual to the boundary magnetic field 
B = B xy {x,r = oo) and the boundary gauge potential Ay{x) = Ay(x,r = oo) of the U(l)v subgroup of 
the SU(Nf) global flavor group of the N = 2 supersymmetric field theory which couples to the boundary 
vector current Jy(x), see, for example, Ref. [37]. For our case, F mrl = 0 but = B xy . 

Rn Eq. 3.1 of [38], T is written in terms of and eu- Here, we have used Eq. 2.35 and 3.4 of [38] 
(which relates th = and = B, respectively) in order to write T explicitly in terms of A h, and 

B. 
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Figure 2 . Critical temperature T C (B) of flavored AT = 4 SYM on R 3 x S ' 1 (3.6) using T° = 
0.175 GeV and X h =23x §*-. 

with increasing number of flavors Nf at zero magnetic field B = 0 and chemical potential 
/x = 0. 

We have plotted (3.6) in Fig. 2 which clearly shows that T C (B ) decreases with increasing 
B <C T 2 in agreement with the inverse magnetic catalysis recently found in lattice QCD 
for B < 1 GeV 2 [19]. 

4 Conclusion 

The observation of inverse magnetic catalysis, for B < 1 GeV 2 , in the recent study of 
QCD on the lattice [19] , instead of magnetic catalysis as predicted by most models of QCD 
[1—16], has posed a serious challenge both for the holographic and non-holographic models 
of QCD. Therefore, motivated by gaining new insight into this problem, we have studied 
the T — B phase diagrams of hard-wall AdS/QCD, and holographic duals of flavored and 
unflavored AA=4 SYM on l 3 x S 1 for B T 2 . 

We have found that the T—B phase diagrams of hard-wall AdS/QCD, and holographic 
duals of flavored and unflavored AA=4 SYM on K 3 x S 1 (which are based on the study 
of the confinement-deconfinement phase transition) are consistent with the recent lattice 
QCD result [19] but opposite to the results of the other holographic models of QCD (which 
are based on the study of the chiral-symmetry breaking phase transition) [11-16]. As it 
can be seen in Fig. 1, Fig. 2, and Eq.(3.3), we have found that the critical temperature 
T C (B) decreases with increasing B <T ] in agreement with the inverse magnetic catalysis 
observed in the recent lattice QCD result for B < 1 GeV 2 [19]. 

However, we would like to emphasize the fact that the vanishing of T C (B), at some 
critical magnetic field B c , as it can be seen in Fig. 1, Fig. 2, and Eq.(3.3), might be the 
artifact of our usage of the backreacted black hole solution only for small magnetic field 
B<T 2 and Nf <C N c limits. Had we used the exact solution for the backreacted black hole 
solution, for any value of B and Nf, the critical temperature T C (B) might not necessarily 
vanish at some critical magnetic field B c . 
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For example, in the large magnetic field B T 2 limit, the Hawking temperature of 
the black hole geometry of [50] is given by T(B > T 2 ) = (A.3), which according to 

the prescription of finding the critical temperature T c given in subsection 3.3, the critical 
temperature T c of the unflavored M = 4 SYM on M 3 X S 1 for large magnetic field B ^$> T 2 
will be T C (B S> T 2 ) = 3 = |T C (H = 0), which means that T C (B T 2 ) > T C (B = 

0). Therefore, T C (B), of unflavored M = 4 SYM on M 3 x S 1 , is not a monotonically 
decreasing function of B, and it will eventually start to increase for large magnetic field B 
and asymptotes to the constant value T C (B T 2 ) = |T C (L> = 0). 

Similarly, for hard-wall AdS/QCD, the critical temperature T C (B 3 > T 2 ) = 

(A.16), but from (3.1), one can see that T C (B = 0) = 2 1 / 4 ^, see also [48, 49], which 
means that T C [B T 2 ) > T C (B = 0). Therefore, T C (B ) of hard-wall AdS/QCD, is not a 
monotonically decreasing function of B , and it will eventually start to increase for large 
magnetic field B and asymptotes to the constant value T C (B T 2 ) = We leave the 

holographic computation of T C (B ) for arbitrary value of B for future investigations. 

This eventual increasing (after initially decreasing for small B ) of T C (B ) for large B 
has also been predicted in some very recent non-holographic studies of the chiral-symmetry 
breaking phase transition in QCD [54, 55], and it has been related to the asymptotic 
freedom and dimensional reduction of QCD at large magnetic field B. 

The inverse magnetic catalysis we found in this paper, from our study of the confinement- 
deconfinement phase transition in the holographic models of QCD (hard-wall AdS/QCD, 
and holographic duals of flavored and unflavored A7=4 SYM on M 3 x S 1 ), is also quali¬ 
tatively similar to the inverse magnetic catalysis observed in the non-holographic studies 
of the confinement-deconfinement phase transition in QCD, such as, MIT bag model [56], 
chiral perturbation theory with two quark flavors [57], and QCD at large- N c [58]. 
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A Confinement-deconfinement phase transition in hard-wall AdS/QCD 

for B > T 2 

In this appendix, we compute the critical temperature T C (B ) of the confinement-deconfinement 
phase transition in hard-wall AdS/QCD for large magnetic field B T 2 . 

For B T 2 , the black hole solution is AdS^ x T 2 black hole or BTZ x T 2 [50] 

ds lh = 3r2 (- f( r )dt 2 + dz 2 ) + -^L ( dx 2 + dy 2 ) + 2 > (A.l) 

M 

f(r) = 1 - ^ . (A.2) 
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and, the Hawking temperature T becomes 

T-i-D'W-fa, (A.3) 

where M = is the mass of the black hole, U{r) = r 2 /(r), the radius of the horizon is 
defined by requiring f(r = 77 J = 0, and (3 is the length of the thermal circle. 

Note that the transverse x — y plane to the direction of the magnetic field B = B z of 
the black hole solution (A.l) is compactified to a 2-torus T 2 [50] which is similar to the 
dimensional reduction that appears in field theories in an external magnetic field B. 

For B 3 > A 2 R ~ rg, the thermal- AdS^ x T 2 solution can be found from the black hole 
solution by setting the mass of the black hole M to zero, that is, 

^thermal = 3r 2 ( — di 2 + dz 2 ) + — j= ( dx 2 + dy 2 ) + 7^77 . (A.4) 

The on-shell Euclidean action Se = Sbh (2.8) for the black hole solution with B ^$> T 2 
(A.l) is 


Sbh — 'S'bulk T Sbndy, 


(A.5) 


where the on-shell Euclidean bulk action S bulk of the black hole solution for B T 2 is 


S'buik = j dt E J dryfg ^4 + -B 2 g xx g yy ^j , 

rP r r \ 

= 6 dt E dr^/g , 

JO Jr h 


(A.6) 

(A.7) 


and, it turned out, we will not need the on-shell Euclidean boundary action of the black 
hole to remove the UV divergences, so we set Sb n d y = 0 for B ^$> T 2 . Hence, the on-shell 
Euclidean action of the black hole Sbh (A.5) becomes 




Sbh = Sbuik + Sbndy = 6 / dt E / dry/g , 


(A.8) 


which (after using the black hole metric for B ^$> T 2 (A.l), evaluating the integrals, and 
simplifying) gives 


S b h = 3 f3B (r\ - r 2 ) 


(A.9) 


The on-shell Euclidean action Se = Sthermai (2-8) for the thermal -AdS solution with 
B 3> Aj R ~ rg (A.4) is 


Sthermai — St bulk T Stbndyj 


(A.10) 


where the on-shell Euclidean bulk action S t b u ik of the thermal-AdS for B A 2 rR ~ rg 


is 


rP' r r a / 2 \ 

Stbuik = J dt E J dry/g ^4 + -^B 2 g xx g yy j , 


(A.ll) 
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and, it turned out, we will not need the on-shell Euclidean boundary action Stbndy of the 
thermal-AdS 1 for B S> A. 2 IR ~ t-q to remove the UV divergences, so we set Stbndy = 0 for 
B ^ ~ r(j. 

The thermal -AdS bulk action Stbuik (A. 11) (after using the thermal -AdS metric for 
B ^ A J R Tq (A.4), evaluating the integrals, simplifying, and using f}' = fiy/J = /3(1 — 

1 r 2 

2 r 2 ) near the boundary) becomes 

Stbuik = 3{i\f]B(r\ - r o) 

= 3/3 B{r\-K 2 h -rl). (A.12) 

Hence, 

^thermal = Stbuik = 3/3 B{r\ - Kl - r%). (A. 13) 

Therefore, ASe (which is the difference between the black hole (A.9) and thermal-AdS’s 
(A. 13) on-shell Euclidean actions) becomes 

A<Sg = <Sbh ^thermal 

= 3 13B (V 2 - r -pj 

= 3(3B (V 2 - ^tt 2 T 2 ) , (A.14) 

for B > T 2 . 

For B 3> T 2 , requiring A S E (T C ) = 0 (A. 14), we find the constraint equation for the 
critical temperature T c to be 

r 2 - ^ 2 T 2 = 0, (A.15) 

which can be solved for T c to give 

T c (B,r o) = ^=-. (A.16) 

\/2 vr 
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